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Granular materials are widely used in industrial processes despite their complex and poorly understood 
mechanical behaviour both in static and dynamic regimes. A prototypical example is the settlement and 
compaction of a granular bed under vibrational loading. The elementary mechanisms of this process are still 
unclear and there is presently no established theory or methodology to predict the settlement and its 
statistical variability. By means of a parametric study, carried out on a full-scale track, and a critical analysis 
of density relaxation laws, we introduce a novel settlement model in coarse granular materials under cyclic 
loading. Our extensive experimental data indicate that the settlement process is governed by three 
independent parameters strongly correlated with the vibration intensity and initial packing fraction. We 
show that the mean settlement is well predicted by the model with its parameter values extracted from 
experimental data. 

Granular materials are both pressure-dependent and density- dependent materials and exhibit a broad 
range of intricate behaviours due to their discrete nature, dissipative interactions and generic 
structural disorder 1 . The packing fraction may vary as a result of particle rearrangements induced 
by shearing or vibrations and it leads to dramatic changes in the structure and mechanical response of a 
granular material 2 " 8 . A long-time logarithmic relaxation law of the packing fraction is systematically observed 
in experiments 910 . In simple compaction models, this behaviour is attributed to the exponentially increasing time 
for the particles to reach a new configuration of lower packing fraction. The case of settlement under cyclic 
loading has, however, been much less investigated. The settlement of granular bed occurs due to both compaction 
and side-wise spreading. An important industrial example is the railway ballast, which undergoes gradual 
settlement under the static and dynamic overloads induced by train traffics 11 " 14 . The readjustment of differential 
settlements requires costly operations on fast-train railways. For this reason, an improved understanding of the 
parameters governing the settlement process is a critical technological challenge for new developments in this 
field. 

In this paper, we show that the total settlement t n under vertical cyclic loading is governed by a logarithmic 
relaxation law as a function of the number N of cycles: 



T N = T a 



1+Bln 1 + 



N 



N 0 J/ 



(1) 



where the three fitting parameters Too, B and N 0 can be evaluated from the loading parameters, namely the 
frequency w (related to the train speed for ballast) and initial packing fraction of the material. Our experimental 
correlations between model and loading parameters show consistently that Too and B depend on the dimensionless 
loading intensity T = (Aco 2 )/(pd 2 /m + g), where A is the vibration amplitude, p is the confining pressure (under 
the sleeper for ballast), d is the average particle diameter, m is the average particle mass and g is gravitational 
acceleration. We also find that the parameter N 0 is linked to the initial packing fraction of the material. In fact, this 
parameter controls the initial settlement rate, and it was systematically determined by means of a light penet- 
rometer in our experiments on ballast material. 
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Figure 1 | Full-scale track model used in the present work. 



Results 

Loading system. We carried out several experimental tests on a full- 
scale track model composed of a granular bed of ballast particles 
(Fig. 1). The bed was subjected to a sinusoidally varying vertical 
overload up to 10 4 cycles for a broad range of amplitudes and 
frequencies. In the same way, the initial state of the granular bed 
was varied and a light dynamic penetrometer was used to 
characterize its packing fraction before each loading test 15 " 17 . This 
device provides the cone penetration strength, which is a parameter 
strongly dependent on the packing fraction of a granular sample. The 
principle of the light penetrometer is to push into the material, 
by manual beating with a hammer, a conic tip of cross-section 
A tip =2 cm 2 . For each blow, the depth of penetration and the input 
energy are measured to calculate the dynamic cone strength q d with 
the corresponding depth using the so-called Dutch formula: 



qd~- 



M 



A tip x eM + P 



(2) 



where M is the mass of the hammer, P is the weight of the driven parts 
during impact (impact head, road and tip), E K is the kinetic energy 
transmitted to the system during the impact and e is the penetration 
per blow. A detailed description of the model and testing apparatus 
can be found in the Methods section below. 

The cumulative settlement may qualitatively be described as 
occurring in three stages as displayed in Fig. 2a. The first stage is 
the initial fast compaction of the granular bed during the first 200 
cycles. Its quasi-linear evolution reflects the uneven distribution of 
voids in the material. The largest voids disappear as a result of local 
particle rearrangements, leading thus to the overall increase of pack- 
ing fraction and settlement. The second stage lasts from 200 to 6000 
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Figure 2 | 
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(a) Settlement of ballast material as a function of the number of cycles with its three stages; (b) Comparison between the experimental 
data obtained from a cyclic loading test with frequency of 3.3 Hz and applied load of 194 X 10 3 N, and several ballast settlement 
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cycles with a nonlinear increase of settlement. This "intermediate" 
stage involves collective particle rearrangements since large voids are 
exhausted in the first stage. Hence, large-scale particle movements 
and longer excursions are necessary for increasingly smaller available 
void space 1,9 . This slow relaxation is quite similar to that observed in 
glasses, spin glasses and flux lattices 18 " 20 . The last stage of settlement 
is the long-term behaviour observed beyond 6000 cycles. In this part 
of the curve, the settlement has a slow linear evolution with the 
number of cycles. No collective particle rearrangements are observed 
so that the slow evolution should be attributed to "rare events" such 
as sudden micro-sliding induced by erosion or grain crushing. 
Beyond this schematic description, we are interested here in a func- 
tional form capable of predicting the settlement for a high number of 
cycles with parameters physically related to and measurable from 
packing fraction and loading conditions. 

Best-fit model. We tested different models proposed for the railway 
ballast settlement based on loading parameters such as vertical stress, 
subgrade stiffness and calibration constants 21 . None of these models 
was found to describe satisfactorily our settlement curves shown in 
Fig. 2b. Shenton expression fits well the beginning of the settlement 
curve but overestimates long-term settlement. The expressions of 
Sato and Hettler do not provide a better fit until 5000 cycles. 
Moreover, a fundamental short- coming of such phenomenological 
models is that they introduce calibration constants that have no 
obvious physical interpretation and cannot be estimated by means 
of independent experimental measurements. 

An alternative approach is the density relaxation law for granular 
materials tested on model systems. This law, called "Chicago fit", was 
obtained by experimental analysis of the evolution of packing frac- 
tion for mono-disperse spherical particles (glass beads) in a cylin- 
drical tube under a series of external excitations consisting of vertical 
shakes or "taps" applied to the container 2,3,9 . This relaxation law is an 
inverse-logarithmic law: 

s.\ Pf-Po 
p(t) = P f-l + Bln(l + t/t 0 ) (3) 



where pf is the ultimate value of packing fraction, p 0 is the initial 
packing fraction, B is a fitting parameter and t 0 is a characteristic 
time. Equation (1) given in Introduction is nothing but a transcrip- 
tion of equation (3) in which the number of cycles AT replaces time, t n 
replaces the packing fraction, and the characteristic time is replaced 
by a characteristic number of cycles N 0 . From 360 independent set- 
tlement data, the parameters for the two models were identified from 
the loading parameters (amplitude and frequency) of each test. 

Table 1 shows the inter- correlations coefficients among the model 
and loading parameters. The parameters Too and B are clearly linked 
to the dimensionless acceleration T whereas N 0 depends on the initial 
packing fraction of the granular material. The correlation between Too 
normalized by the initial thickness of the granular layer H 0 and T has 
a high correlation coefficient R 2 = 0.9 (Fig. 3a), thus suggesting a 
linear relationship between the two parameters. Different values of 
amplitude and confining pressure lead to a similar behaviour with 
the predicted values of tJH 0 . In the same way, for different values of 
A, co andp but similar values of T, the settlement curves show a quite 
similar behaviour, as seen in Fig. 4a. The parameters A, co are inde- 
pendent from each other, while T increase as a function of both 
(Fig. 4b and 4c). Y appears thus to be the governing parameter of 
settlement. 

Figure 3b shows the correlation between B and T. Despite the 
relative dispersion of data due to the inherent variability of the sys- 
tem, it is seen that B globally declines as the acceleration parameter 
increases. Regarding the N 0 parameter, a logarithmic relationship 
with the tip strength (normalized by p) is observed; see Fig. 3c. 
Indeed, increasing the initial packing fraction, in correlation with 
the penetration strength q d , leads to higher values of N 0 and therefore 
lower initial slope or slower settlement in the first loading cycles. 

In Fig. 5a the experimental data are compared with the above 
model of settlement. This model describes quite well the evolution 
of settlement up to a high number of cycles and for the three stages of 
settlement when the three model parameters are adjusted. In order to 
assess the quality of this model, we calculate the mean square error 
(MSE) obtained between all the experimental and predicted values of 
settlement. MSE is defined as 

MSE(6\6)=((6-6f^ (4) 

where 0 is the vector of predicted values for each settlement test and 0 
is the vector of the corresponding measured values. In Fig. 5b we have 
plotted the cumulative distribution function of MSE between 360 
experimental tests: 40% of the settlement curves are predicted with 
an error below 5% and in 60% of tests the error is below 10%. This 
indicates that the model based on the Chicago fit provides a quite 
accurate estimation of the evolution of settlement in coarse granular 
materials and, what is more, the model parameters have a clear 




0.05 



1.6 
I. 
1. 
L 



-0.8 
0.6 
0.4 
0.2 

°(fc>0 



• • data 
■-- model 



•-..V 



310 



0.01 



0.02 0.03 

r 



0.04 0.05 




3.5 4.0 



Figure 3 | Correlations between the model parameters and loading parameters, (a) The relation between settlement normalized by the thickness of the 
granular bed xJH 0 and T; (b) The relation between B and T; (c) The relation between N 0 and qjp. 
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(a) Settlement for different values of A, co andp, but the same value of T; (b) Values of T as a function of the amplitude for co= 18.85 Hz and 
10 3 Pa; (c) Values of T as a function of co, for A= 0.8 X 10" 3 m and p= 170 X 10 3 Pa. 



correlation with the packing fraction and loading parameters and can 
be obtained by independent tests. 

Discussion 

In this work, we carried out a parametric study on a full-scale ballast 
track in order to test various models for the prediction of settlement 
in coarse granular materials under cyclic forcing. The settlement as a 
function of the number of cycles occurs in three stages with nearly the 
same functional dependence but varying with the initial state of the 
material and vibration amplitude. We showed that a model based on 
the Chicago density relaxation law provides an accurate estimation of 
settlement for a high number of cycles with parameters that reflect 
the loading parameters and initial mechanical state of the material 
and can be measured by independent tests. This model, together with 
its testing procedure, provides a reliable approach in applications 
where long-term settlement needs to be estimated in coarse granular 
materials. 

Our findings raise also a fundamental question as to the predict- 
ability of the mechanical response of granular materials despite their 
natural variability and complex shapes of ballast grains. The beha- 
viour is chaotic in the sense that small variations in the initial con- 
figuration of a granular packing are expected to be amplified with the 
number of loading cycles. This is not what we observe although 
fluctuations are observed between independent tests as shown in 
Fig. 3. This indicates that the fluctuating parts of consecutive 



incremental settlements average out and the system tends to a 
well-defined ultimate mechanical state. In other words, the mech- 
anical state of the granular material under vibrations is determined 
by its distance from the ultimate state. This is consistent with the 
form of equation (1) where the relevant settlement parameter 
appears to be t^/Toc with Too only depending on the dimensionless 
acceleration parameter T. This feature of compaction can be generic 
to all systems characterized by jamming transition and in which a 
maximally jammed state may be defined 19 ' 20 . Paradoxically, in such 
systems the long-term behaviour may be more accurately predicted 
than the short-term behaviour. 

Methods 

The track model is composed of two monoblock sleepers lying on a ballast layer of 
thickness 0.35 m. The ballast layer is confined in a wood box of dimensions 4.1 1 m by 
1.5 m. In this track model, the subgrade is represented by a wooden plate. The 
sleepers are fixed together by a rail section. A hydraulic cylinder applies the sinusoidal 
loading on the track. With the aim of providing a uniform load distribution, a beam is 
fixed to a hydraulic cylinder applying the overload on the rails. To study the influence 
of loading on settlement, we apply a sinusoidal load during 10 4 cycles. The peak load is 
varied from 194 to 272 X 10 3 N and the frequency from 3.3 Hz to 6 Hz. The cyclic 
loading is applied by a hydraulic cylinder on two sleepers placed on the granular bed. 
Before applying the cyclic loading, we characterized with several light penetration 
tests the initial mechanical state of the granular bed in our full-scale ballast track. In 
the course of cyclic loading, we measured with 4 non-contact optical sensors the 
evolution of settlement as a function of the number of cycles. The tests were carried 
out with 15 different combinations of the mean load and frequency and 6 inde- 
pendent initial conditions. 
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Figure 5 | (a) Prediction of settlement by a logarithmic law in comparison to the experimental data for two tests (frequency of 4.5 Hz and applied load of 
239 X 10 3 N; Frequency of 3.3 Hz and applied load of 194 X 10 3 N). The inset shows the data on a log-linear scale; (b) Cumulative distribution function 
of the mean square error (MSE) of model prediction compared to 360 experimental tests. 
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